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INTERNATIONAL  CONFERENCE  ON  INDUSTRIAL 
AND  APPLIED  MATHEMATICS 


l  INTRODUCTION 


Four  international  societies  of  ap¬ 
plied  mathematics — die  Gesellschaft  fiir 
Angewandte  Mathematik  and  Mechanik 
(GAMM) — of  the  German-speaking  countries, 
the  Institute  of  Mathematics  and  its  Ap¬ 
plications  (IMA)  of  the  UK,  the  Society 
for  Industrial  and  Applied  Mathematics 
(SIAM)  of  the  US,  and  Soci£t£  de  Math£- 
matiques  Appliqu£es  et  Industr ielles 
(SMAI)  France,  joined  forces  in  Paris 
from  June  29  through  July  3  to  conduct 
the  First  International  Conference  on 
Industrial  and  Applied  Mathematics. 

The  popularity  of  the  meeting  ex¬ 
ceeded  all  expectation,  with  more  than 
1500  attendees,  more  than  1000  papers 
(which  were  given  in  17  parallel  ses¬ 
sions),  and  130  poster  sessions.  The 
topics  ranged  from  vorticltv,  chaos,  dy¬ 
namical  systems,  robotics,  and  control 
to  optimization,  parallel  processing, 
multigrid  methods,  algorithm  development, 
stochastic  modeling,  lattices,  and  net¬ 
working.  In  addition  to  the  parallel  ses¬ 
sions,  there  were  plenary  sessions  which 
featured  a  total  of  16  invited  speakers. 

The  country  with  the  greatest  number 
of  attendees  (as  might  be  expected)  was 
France,  with  310,  followed  by  the  US 
(.’r*0)  ,  West  Germany  (160),  and  UK  (110). 
The  remainder  of  the  attendees  came  from 
45  other  countries  in  North,  Central, 
and  South  America,  Africa,  Western  and 
Eastern  Europe,  the  USSR,  Asia,  and 
Australia. 

I  will  be  reporting  on  papers  given 
in  the  area  of  fluid  mechanics. 


L  VORTEX  DOMINATED  FLOWS 


In  many  high  Reynolds  high  number 
problems,  vortices  play  a  dominant  role 
in  establishing  the  behavior  of  the  flow. 
Examples  are  the  role  played  by  the  lead¬ 
ing  edge  vortex  in  producing  a  lift  aug¬ 
mentation  on  delta-wing  aircraft  and  the 


role  of  vortices  in  explaining  and  con¬ 
trolling  the  rate  of  mixing.  Two  con¬ 
trasting  approaches  are  possible  in  r.r- 
merical  vortex  simulations.  One  is  the 
so-called  Vortex  Dynamics  approach  in 
which  the  vorticies  are  treated  discrete¬ 
ly  in  the  context  of  a  vorti.cJ.ty/ stream 
function  formulation  of  the  incompres¬ 
sible  Navier-Stokes  equations  and  the 
vortex  paths  are  traced  by  particle 
tracking  methods.  The  other  (less  mathe¬ 
matically  sophisticated)  approach  is  to 
"capture"  the  vorticies  in  the  course 
of  more  or  less  conventional  solutions 
of  the  incompressible  or  compressible 
Navier-Stokes  or  Euler  equations  in  which 
sufficient  local  grid  refinement  has  beem 
introduced  to  resolve  the  vortr::  pat¬ 
terns  . 

In  an  invited  paper,  A.  Mnjda  cf 
Princeton  University  (New  Jersey)  re¬ 
viewed  the  applications  of  the  vorfex 
dynamics  approach  and  presented  some  verv 
recent  results  involving  large-scale  com¬ 
putations  using  this  method.  Specifical¬ 
ly,  he  reported  on  calculations  which  he 
had  done  to  model  the  evolution  of  the 
primary  vortex  on  an  elliptical lv  loaded 
delta  wing  (with  and  without  the  effects 
of  the  fuselage)  and  to  model  the  rola  up 
of  the  KIrchoff  vortex  in  a  planar  vortex 
sheet  subjected  to  a  periodic  perturba¬ 
tion.  Embodied  in  these  calculations  are 
what  he  termed  "real  breakthroughs"  which 
involved  (I)  the  use  of  streamlines  o: 
the  Inviscid  theory  to  obtain  nearly  op¬ 
timal  convergence  properties  and  il’)  the 
use  of  multiple  expansions  to  reduce'  tin- 
operator  count  in  the  vortex  summation 
operation  to  order  N,  where  N  is  the  num¬ 
ber  of  vorticies  used  in  the  simulation. 

Incredibly  complicated  vortex  pat¬ 
terns  are  produced  when  vorticies  of  dif¬ 
ferent  sign  are  present  such  as  in  a  del¬ 
ta  wing  calculation.  This  is  illustrated 
in  Figure  1,  which  shows  the  cross-flow 
velocity  contours  of  a  delta  wing  at  a- 
angle  of  attack  of  10°.  The  calculation 
was  done  for  a  free  stream  Mach  number  of 
2  with  the  additional  complication  that 
the  wing  was  in  rolling  oscillation.  The 
results  were  obtained  by  O.A.  Kandil  (Old 
Dominion  University,  Norfolk,  Virginia) 
not  with  a  vortex  dynamics  method  but  bv 
means  of  a  time-accurate  solution  of  the 
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Euler  equations.  These  calculations  are 
particularly  interesting  since  it  is  the 
::rst  time  that  such  calculations  have 
been  attempted  for  a  maneuvering  wing. 
The  calculations  feature  a  three-dimen¬ 
sional,  fully-vectorized  computer  program 
which  solves  the  Euler  equations  for  the 
relative  motion  in  a  moving  frame  of  ref¬ 
erence  using  central-difference,  finite- 
volume  discretization,  and  a  four-stage 
Eunge-Kut ta  integration  scheme.  Several 
possibilities  exist  lor  t  he  choice  of 
variables  and  (in  the  ease  ol  subsonic 
I  lew)  the  manner  in  which  the  condition 
of  constant  entropy  is  introduced. 


Kandll  found  that  only  the  algebraic 
entropy  condition  gave  the  proper  (un¬ 
separated)  leading  edge  solution  for  a 
rounded  edge  delta  wing  while  for  sharp- 
edged  delta  wings  it  didn't  make  any  dif¬ 
ference  whether  the  constant  entropy  con¬ 
dition  was  introduced  algebraically  or  a 
constant  value  was  prescribed  by  means  of 
a  partial  differential  equation  of  en¬ 
tropy  transport.  In  solving  a  similar 
problem  (but  tor  a  nonmaneuvering  case) 
C.-ll.  Hriineau  of  the  Uuivcrslte  de  I’.irls- 
Sud  (Kranic)  I  ound  that  in  order  to  get  .1 
stable  solution  ol  the  steady  Euler  equa¬ 
tions  for  both  subsonic  and  supersonic 


flows,  he  needed  to  solve  an  explicit 
energy  transport  equation.  If,  instead, 
a  constant  value  of  the  stagnation  en¬ 
tropy  was  enforced  by  means  of  an  alge¬ 
braic  relationship,  convergence  could 
only  be  obtained  for  subsonic  free-stream 
t  low. 

By  far  the  most  extensive  vortex 
wing  computations  were  reported  by  A. 
Kizzi  of  The  Aeronautical  Research  Insti¬ 
tute  of  Sweden.  He  carried  out  both  Euler 
and  Navier-Stokes  calculations  using  his 
well-known  finite-volume,  cell-centered 
scheme  with  one-step,  three-stage  Runge- 
Kutta  Integration.  He  characterized  his 
Euler  calculations  as  being  disappoint¬ 
ing.  He  expected  from  the  outset,  of 
course,  that  the  Euler  calculations  could 
not  capture  the  secondary  separation  re¬ 
sulting  from  the  cross-flow  shock.  How¬ 
ever,  it  was  expected  that  the  overall 
character  of  the  solution  as  evidenced 
by,  say,  the  static  pressure  distribution 
on  the  wing  would  be  more  or  less  cor¬ 
rect.  This  turned  out  not  to  be  the  case. 
The  static  pressure  rise  due  to  the  pri¬ 
mary  vortex  seemed  to  be  too  strong  and 
located  too  far  outboard  on  the  wing,  a 
situation  which  did  not  improve  even 
when  more  than  1  million  points  were  used 
in  the  calculation.  Additional  details 
of  Rizzi 's  calculations  can  he  found  in 

41-1:22  [1987]). 

Rizzi  has  just  obtained  the  first 
results  from  his  new  Navier-Stokes  code, 
which  uses  the  same  cell-centered  scheme 
as  the  inviscid  calculations.  The  deriva¬ 
tives  appearing  in  the  viscous  terms  were 
evaluated  by  means  of  a  single-grid  two- 
'-tep  application  of  the  gradient  theorem 
CPevret  and  Taylor,  1983).  A  Mach  0.85 
case  at  10°  Incidence  was  run  as  a  lami¬ 
nar  calculation  at  a  Reynolds  number  of 
2.4<10f'  employing  400,000  cells,  100,000 
of  which  were  In  the  boundary  layer.  Sev¬ 
eral  qualitative  differences  from  the 
earlier  inviscid  calculations  were  im¬ 
mediately  apparent: 

o  The  origin  of  the  primarv  vortex  was 
at  the  wing  apex  rather  than  along  the 
leading  edge. 

o  The  Navier-Stokes  equations  picked  tip 
the  secondary  separation  and  perhaps 
even  a  tertiary  one. 


o  A  cascade  of  vorticies  seemed  to  arise 
in  the  vicinity  of  the  (rounded"' 
leading  edge  of  the  wing. 

The  most  striking  difference,  however, 
was  in  the  greatly  improved  accuracy  ot 
the  predicted  static  pressure.  The  new 
origin  of  the  primary  vortex  resulted  in 
the  pressure  peaks  being  moved  inboarc 
the  proper  location,  and  the  level  of  the 
pressure  rise  was  significantly  (50  per¬ 
cent)  less  than  the  inviscid  calculation 
and  in  excellent  agreement  with  the  ex¬ 
periment.  Rizzi  was  pleased  with  these 
improved  results  and  also  with  the  fact 
that  the  addition  of  the  viscous  terras 
increased  the  CPU  time  by  only  50  per¬ 
cent.  An  important  conclusion  which  can 
be  drawn  from  this  work  is  that  viscous 
effects  play  a  much  more  important  role 
in  establishing  the  proper  qualitative 
and  quantitative  behavior  of  delta  wing 
flows  than  was  previously  thought  and,  in 
lact,  are  absolutely  essential  (at  least 
for  rounded  leading  edge  wings)  in  order 
to  establish  the  proper  surface  flow 
topology  and  static  pressure  riistribu- 
t  ion. 

Rizzi  then  turned  his  attention  to 
the  detailed  structure  of  the  primary 
vortex,  as  revealed  by  his  high-resolu- 
t  ion  (more  than  1  million  point)  Fuler 
calculation.  The  fundamental  problem 
which  arises  in  attempting  to  extract  the 
detailed  structure  of  the  vortex  from  an 
huler  calculation  is  that,  as  the  core  of 
the  vortex  is  approached,  it  becomes  im¬ 
possible  to  control  the  amount  of  numer¬ 
ical  error  which  is  present  (due  to  the 
velocity  gradient,  which  approaches  in¬ 
finity  at  the  core  of  the  vortex).  Thus, 
a  truly  inviscid  simulation  of  the  vortex 
structure  is  impossible.  Near  the  outer 
boundary  of  the  vortex,  however,  the 
Euler  calct  lations  are  relatively  tree 
from  numerical  dissipation  and  a  more  or 
less  clearly  defined  vortex  sheet  can  *^e 
distinguished  from  the  calculations.  As 
rhe  vortex  sheet  spirals  inward,  however, 
the  Identity  of  the  sheet  is  ’ost  and 
diffused  core  is  produced,  the  calcula¬ 
tion  of  which  is  dominated  bv  the  effects 
>>f  numerical  viscosity. 

It  is  surprising  that  the  less 
in  total  pressure  predicted  bv  Rizzi V 


i 


calculations  agrees  so  well  with  the  ex¬ 
perimentally  determined  value.  In  a  cal¬ 
culation  of  a  75°  swept  wing  at  an  incid¬ 
ence  of  10°  and  a  Mach  number  of  1.05, 
Rizzi  found  a  total  pressure  loss  of  70 
percent  compared  with  the  experimental 
value  of  75  percent.  Moreover,  the  mag¬ 
nitude  of  the  loss  in  total  pressure  ap¬ 
peared  to  be  independent  of  the  level 
',ic  numerical  dissipation,  a  result  subr.e- 
ouentlv  confirmed  bv  Murraan  (FSIf  41-1:22 
r 1 987 ’ )  . 

3  LARGE-SCALE  COMPUTING  OF  AEROSPACE 
P  ROPI.EMS 

In  separate  but  reflated  presenta¬ 
tions,  P.  Perrier  and  J.Periaux,  both  of 
Avions  Marcel  Dassault-Br^guet  Aviation 
w\MD-!'A),  reviewed  recent  progress  in 
algorithm  development  and  computer  design 
which  have  allowed  computational  fluid 
dynamics  (CFP)  to  assume  the  vital  role 
which  it  now  plays  in  the  design  of  to¬ 
day's  aerospace  vehicles.  In  20  years, 
FD  has  grown  from  zero  to  nearly  50  per¬ 
cent  of  the  development  cost  of  new  aero¬ 
space  vehicles.  We  have  now  reached  the 
point  where  10 operations  can  be  com¬ 
pleted  In  a  10-hour  day  and  this,  he 
said,  is  the  threshhold  at  which  we  can 
begin  to  obtain  accurate  simulations  of 
the  flow  fields  about  physically  realis¬ 
tic  aerospace  vehicles.  Thanks  to  the 
strong  interaction  which  has  existed  at 
AVD-BA  with  fundamental  research  organi¬ 
zations,  such  as  Office  National  d' Etudes 
et  de  Kecherches  Alrospatiales  (ONERA) , 
Institut  National  de  Recherche  en  In- 
formatique  et  en  Automatique  (INRIA),  and 
various  European  universities,  fast  and 
accurate  solution  of  large  sets  of  non¬ 
linear  partial  differential  equations  Is 
now  possible.  Although  the  details  of 
the  technique  vary  depending  upon  the 
application  (aircraft,  turboaachlnery , 
missiles,  or  reentry  vehicles),  it  is 
AMP-BA's  position  that  the  only  method 
which  Ik  capable  of  Che  drastic  mesh 
adapt  ions  required  to  economically  simu¬ 
late  relatistic  aerospace  vehicle  flow 
fields  Is  a  i. on-structured  finite  element 
algorithm  using  domain  decomposition  and 
multigrid  acceleration. 


Periaux  described  several  impressive 
examples  of  aerospace  applications  Illus¬ 
trating  the  current  status  of  CFD  tech¬ 
nology  beginning  with  the  calculation  of 
potential  flow  over  a  complete  transport 
aircraft  done  more  than  10  years  ago, 
two-  and  three-dimensional  Navier-St okes 
simulations  ol  transfer  inlet  flow  calcu¬ 
lations  of  the  flow  over  a  complete  mili¬ 
tary  and  commercial  aircraft,  and 
lations  of  three-dimen?;  f  ona  1  reentry 
bodies.  Despite  the  convincing  nature  ol 
the  calculation  of  the  local  static  pres¬ 
sure,  for  example,  upon  closer  examina¬ 
tion  the  prediction  of  other  quantities 
is  not  as  satisfactory.  For  example,  he 
showed  the  calculation  of  the  flow  in  the 
vicinity  of  the  wing  root  of  a  subsonic 
transport  which  exhibited  an  unwanted, 
and  physically  unrealistic,  entropy  in¬ 
crease.  This  can  be  traced  to  inadequate 
grid  resolution.  Similar  difficulties 
occur  in  attempts  to  model  the  curvature, 
of  engine  inlets  and  exhaust  nozzles, 
Periaux  concluded  that  although  there  is 
good  reason  to  be  excited  about  the  de¬ 
sign  of  possibilities  of  today's  CFD 
methods,  much  research  remains  to  be  done 
on  development  of  accurate  algorithms  to 
make  optimum  use  of  t he  current  genera¬ 
tion  of  supercomputers. 

In  the  US,  finite  volume  codes  are 
much  more  common  than  finite  element 
codes.  S.  Chakravarthv  of  Rockwell  In¬ 
ternational  (US)  presented  the  results 
of  his  Euler  and  Navier-Stokes  calcula¬ 
tions  using  a  class  of  numerical  methods 
based  on  total  variation  diminishing 
(TVD)  formulation.  In  the  past  12  months 
he  has  Imbedded  this  method  into  a  uni¬ 
fied  code  which  has  the  capability  of 
doing  multizone,  multigrid,  and  either 
explicit  or  implicit,  time  or  space 
marching  calculations.  The  results  of  his 
three-dimensional  Euler  calculations  for 
supersonic  flow  over  a  complete  shuttle/ 
external  tank/solid  rocket  booster  assem¬ 
bly  produced  static  pressures  which 
agreed  very  well  with  flight  test  data. 
Calculations  were  also  shown  for  a  gener¬ 
ic  hypersonic  transport  configuration 
which  involved  the  solution  of  the  two- 
dimensional  Navier-Stokes  equations.  The 
calculations  were  done  for  a  Mach  number 
of  25  and  featured  a  combination  of  time 
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marching  at  the  nose,  patched  to  space 
marching  for  the  remainder  of  the  flow 
field.  Having  now  obtained  a  comprehen¬ 
sive  code  for  aerodynamic  simulation, 
Chakravarthy 's  next  step  is  to  combine  it 
with  propulsion,  structure,  acoustics, 
and  radar  cross  section  codes  in  order  to 
develop  a  truly  comprehensive  CFD  design 
package . 

J.  Chattot  of  MATRA  (France)  de¬ 
scribed  work  undertaken  to  meet  a  much 
more  modest  objective — the  development 
of  an  aerodynamic  design  code  for  high 
angle  of  attack  missile  operation.  The 
basis  of  Chattot's  code  is  an  unsteady 
finite-volume  structural  mesh,  Euler  cal¬ 
culation  structured  meshes  employing  Lax- 
Wendroff-Ni  integration  procedure.  Use 
of  an  explicit  integration  procedure, 
Chattot  said,  was  motivated  by  his 
desire  to  facilitate  the  use  of  multi¬ 
zone  solution  techniques.  As  in  Chakra- 
varthy's  code,  multigrid  was  used  to  en¬ 
hance  convergence  rate.  The  results 
showed  that  normal  force  and  pitching 
moments  were  well  predicted  for  models 
of  considerable  geometric  complexity  up 
to  an  angle  of  attack  of  15°.  To  improve 
the  predictions,  a  collaborative  effort 
is  being  undertaken  with  ONERA  to  couple 
the  code  with  an  unsteady  boundary  layer 
calculation.  Further  details  concerning 
this  work  can  be  found  in  (FSN  41— 
4:215  [1987]). 

ONERA  has  also  developed  a  series  of 
Euler  codes  for  the  simulation  of  complex 
three-dimensional  flows  past  aerodynamic 
configurations  based  on  the  Lax-Wendroff 
procedure.  This  work  was  described  by  P. 
Morice  of  ONERA-Chatillon.  Each  of  the 
codes  is  well  suited  to  a  specific  field 
of  application,  namely,  low  speed,  tran¬ 
sonic,  or  hypersonic  flows.  Application 
of  the  new  three-dimensional  Euler  code 
to  the  calculation  of  a  unperson  ic  ellip¬ 
tical  jet  using,  jet  surface  "fitting"  was 
described.  In  these  calculations  the  sur¬ 
face  of  the  jet  was  moved  in  accordance 
with  appropriate  characteristic  relation¬ 
ships  in  order  to  match  the  prescribed 
constant-pressure  boundary  condition.  A 
surprising  amount  of  three-dimensionality 
was  observed  in  the  calculations  for  the 
relatively  modest  eccentricity  of  the 
nozzle  cross  section. 


For  supersonic  and  hypersonic  appli¬ 
cations  ONERA  has  developed  a  second  or¬ 
der  of  accurate,  explicit,  upwind  (Van 
I,eer-type)  flux-split,  nodal  scheme.  At 
the  present  time,  real-gas  equilibrium 
calculations  are  included,  and  work  is 
underway  to  add  the  effects  of  chemical 
reaction.  Two  applications  were  shown, 
one  to  a  rectangular  supersonic  jet  ex¬ 
hausting  into  a  supersonic  cross  stream 
and  another  involving  a  missile  at  a  45° 
negative  angle  of  attack.  For  the  jet 
calculation,  a  35x40x45  grid  was  used, 
which  revealed  a  considerable  amount  of 
fluid  mechanics  detail  in  the  free- 
stream/jet  interaction  region.  In  partic¬ 
ular,  the  expansion  of  the  free  stream 
around  the  back  of  the  jet  was  well  re¬ 
solved,  as  was  the  cross-stream  flow,  as 
evidenced  by  cross-flow  velocity  contours 
plotted  at  several  downstream  occasions. 
The  grid  refinement  procedure  appeared  to 
be  far  less  sophisticated  than  that  con¬ 
tained  in  AMD-BA's  finite  element,  codes. 
Morice  reported  that  work  is  presently 
underway  to  add  such  a  capability  to  the 
program  as  well  as  to  develop  improved 
turbulence  modeling  and  transition  cri¬ 
teria  for  the  eventual  viscous  extension 
of  these  calculations.  The  objective, 
like  that  described  by  Chakravarthy,  is 
to  develop  reliable  aerodynamic  simula¬ 
tion  codes  for  industrial  design  use. 

4  CONVERGENCE  ACCELERATION  METHODS 

Perhaps  the  most  popular  convergence 
acceleration  method  these  days  is  the 
multigrid  method.  To  set  the  stage  for 
subsequent  sessions  or  multigrid  theory 
and  application,  a  plenary  session  was 
organized  'n  which  an  invited  paper  re¬ 
viewing  recent  progress  in  inultigrid 
theory  was  presented  by  W.  Hackbusch  of 
the  Christian  Albrechts  University  (Kiel, 
West  Germany) . 

The  problem  with  using  a  single 
grid,  Hackbusch  said,  is  that  low  fre¬ 
quency  errors  in  the  numerical  solution 
are  difficult  to  remove.  The  essence  of 
the  multigrid  concept  is  that  an  error  of 
wavelength  A  is  most  easily  eliminated  on 
a  mesh  of  size  h,  where  A  is  approximate¬ 
ly  equal  to  h.  In  view  of  this,  a  giver, 
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flow  problem  Is  solved  on  a  sequence  of 
grids  in  which  the  low-frequency  error, 
when  viewed  from  the  perspective  of  a 
grid  of  increased  coarseness,  appears  to 
be  of  higher  frequency  and  therefore, 
easier  to  remove.  By  subsequent  repeti¬ 
tions  of  a  process  in  which  a  solution  is 
sought  on  coarse  grids  and  the  results 
transferred  back  to  serve  as  the  first 
approximation  for  the  solution  on  the 
next  finest  grid  for  the  next  sequence  of 
calculations,  dramatically  enhanced  con¬ 
vergence  rates  can  be  achieved.  Although 
this  concept  has  been  around  for  some 
time,  it  has  been  implemented  in  a  modern 
computational  context  by  A.  Brandt  and 
placed  on  a  firm  theoretical  basis  only 
within  the  past  decade.  The  past  several 
years  have  seen  the  application  of  multi¬ 
grid  methods  to  the  solution  of  elliptic, 
parabolic,  and  hyperbolic  problems,  sin¬ 
gular  perturbation  problems,  and  integral 
equations  of  the  second  kind.  Of  partic¬ 
ular  interest  to  Hackbusch  is  the  devel¬ 
opment  of  robust  methods  where  conver¬ 
gence  enhancement  effects  are  achieved 
independent  of  mesh  spacing  and  the  geom¬ 
etry  of  the  computational  region.  Par¬ 
ticular  difficulties  occur  in  the  appli¬ 
cation  of  multigrid  methods  to  two-dimen¬ 
sional  regions  where  high  levels  of  ani¬ 
sotropy  occur  in  both  coordinate  direc¬ 
tions.  For  this  he  proposed  the  use  of  a 
new  multigrid  method  and  a  novel  prolon¬ 
gation  method  employing  four  contiguous 
grids . 

A  large  number  of  the  presentations 
focused  on  various  aspects  of  multigrid 
refinements,  such  as  the  relative  effec¬ 
tiveness  of  various  solution  techniques 
(called  "smoothers"  in  multigrid  termi¬ 
nology)  and  various  methods  for  trans¬ 
ferring  coarse  grid  solution  into  the 
next  finer  mesh  ("prolongation").  The 
most  interesting  engineering  application 
of  the  multigrid  method  was  presented  by 
P.  H.  Gaskell  of  the  University  of  Leeds 
^Uk).  He  described  a  multigrid  calcula¬ 
tion  for  the  flow  in  a  lid-driven  and  a 
thermally  driven  cavity  using  a  finite 
volume,  staggered  grid  approach.  In  ad¬ 
dition  to  showing  the  effectiveness  of 
multigrid  (a  time  savings  of  more  than  a 
factor  of  100  was  realized)  the  calcula¬ 
tions  showed  the  superiority  of  the  quad¬ 


ratic  upwind  interpolation  (QUICK)  scheme 
in  contrast  with  the  hybrid  algorithm;  it 
also  showed  the  advantage  of  using  sym¬ 
metrical  coupled  Gauss-Seidel  (SCGS)  it¬ 
erations  in  comparison  with  the  uncoupled 
iteration  technique  in  which  an  independ¬ 
ent  "pressure-correction"  calculation  is 
made.  Gaskell  found  that  multigrid,  in 
combination  with  SCGS  and  the  QUICK  dis¬ 
cretization  scheme,  allowed  him  to  attain 
the  nearly  linear  dependence  of  computa¬ 
tional  speed  on  the  number  of  grid  points 
that  one  would  expect  from  a  classical 
multigrid  method.  Unfortunately,  in  re¬ 
gions  of  high  shear  (for  example,  near 
the  moving  wall  in  the  lid-driven  cavity 
case)  the  unbounded  characteristic  of  the 
QUICK  scheme  causes  difficulty.  This  was 
seen  in  calculations  of  a  turbulent  lid- 
driven  cavity  case  where  unphysical  os¬ 
cillations  produced  negative  values  of 
the  turbulence  kinetic  energy  and  would 
have  caused  the  solution  to  diverge  if  no 
corrective  action  been  taken.  To  overcome 
this  difficulty,  Gaskell  has  developed 
what  he  calls  his  Curvature  Compensated 
Convective  Transport  (CCCT)  discretiza¬ 
tion.  It  fulfills  what  he  has  identified 
as  the  minimum  and  necessary  requirement 
for  boundedness  which  he  has  called  the 
Convection  Boundedness  Criterion  (CBC) . 
Unlike  hybrid  discretization ,  however. 
CCCT  assures  a  bounded  solution  without 
the  overwhelming  numerical  dissipation 
which  is  associated  with  the  hybrid  meth¬ 
od.  Further  information  about  this  work 
can  be  found  in  Gaskell  and  Lau  (1988). 
Currently,  work  is  underway  to  replace 
the  point-by-point  SCGS  solver  with  a 
line-by-line  technique  and  the  use  of 
quadratic  interpolation  for  the  prolonga¬ 
tion  and  restriction  operators  which 
would  further  improve  the  performance  of 
the  technique. 

S  CALCULATION  METHODS  FOR  TURBULENT  FLOWS 


Two  classical  methods  for  calcu¬ 
lating  turbulent  flows  without  the  use 
of  a  turbulence  model  are  direct  and 
large  eddy  simulation.  Of  the  two  meth¬ 
ods,  direct  simulation  is  by  far  the  most 
computationally  intensive  since  it  has  as 
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its  objective  the  numerical  resolution 
(insofar  as  possible)  of  all  relevant 
turbulent  scales.  N.  Satofuka  of  the 
Kyoto  Institute  of  Technology  (Japan) 
described  his  direct  simulation  calcula¬ 
tions  which  employed  the  vorticity/stream 
function  formulation  of  the  Navier-Stokes 
equations.  For  the  discretization  of  the 
spatial  derivatives  in  the  vorticity 
transport  equation  Satofuka  used  the 
method  of  lines.  For  the  solution  of  the 
vorticity  equations  he  used  a  4th-order- 
accurate  Runge-Kutta  scheme,  and  for  the 
solution  of  the  stream  function  equation 
he  used  a  spectral  method.  Satofuka  was 
particularly  interested  in  the  decay  of 
two-dimensional  homogenous  isotropic  tur¬ 
bulence.  In  comparing  the  results  of  his 
totally-explicit  calculations  with  the 
more  familiar  pseudo-spectral  calcula¬ 
tions,  he  found  that  the  vorticity  dis¬ 
tribution  maps  were  identical  and  that 
the  explicit  calculations  required  from 
only  1/4  to  1/5  of  the  time  required  for 
the  pseudo-spectral  calculations .  His 
calculations,  which  were  carried  out  at 
a  Reynolds  number  of  25,000  and  60,000 
were  made  with  meshes  of  128*128,  512*512 
and  1024*i024  points. 

The  second  approach,  where  the  small 
scales  are  modeled  and  only  the  large- 
scale  turbulence  is  calculated,  is  called 
large  eddy  simulation.  An  example  of 
such  a  calculation  was  presented  by  J.P. 
Chollet  of  the  Institut  de  M£canique  de 
Grenoble  (France).  The  first  case  which 
Chollet  described  was  a  two-dimensional 
incompressible  plane  mixing  layer.  To 
visualize  the  flow,  he  used  both  isovor- 
ticity  maps  and  color  contours  illus¬ 
trating  the  convection  of  a  passive  scal¬ 
er.  To  simulate  the  effects  of  three- 
dimensional  disturbances,  which  were,  of 
course,  absent  in  this  two-dimensional 
calculation,  he  introduced  random,  small- 
scale  perturbations  in  the  subgrid  scale 
model.  These  were  responsible  for  initi¬ 
ating  the  Kirchof f-Taylor  instability  in 
the  mixing  layer  from  which  there  fol¬ 
lowed  the  appearance  of  discrete  vortices 
and  subsequent  vortex  pairing.  Develop¬ 
ment  of  an  appropriate  subgrid  scale 
model  is  by  no  means  a  closed  issue  and 
Chollet  presented  two  methods  for  doing 
this,  one  which  made  use  of  a  5/3  power 


decay  of  turbulence  kinetic  energy  with 
wave  number,  and  the  other  (and  far  more 
computationally  intensive)  used  the  EDQNM 
closure  model.  In  his  review  of  the  re¬ 
normalization  group  (RNG)  theorv,  V. 
Yakhot,  Princeton  University  (New  Jersey) 
indicated  his  support  of  the  large  eddy 
simulation  philosophy.  He  said  that  lie 
believed  that  no  one  could  ever  afford 
to  compute  a  direct  simulation  of  an 
aerodynamic  flow  at  full-scale  Reynolds 
number.  RNG,  he  believes,  offers  an 
ideal  subgrid  scale  model  for  large  eddy 
simulations.  RNG  allows  the  evaluation 
of  transport  coefficients  and  transport 
equations  for  the  large-scale  (small  wave 
number)  modes  by  dynamic  scaling  and  in¬ 
variance  together  with  iterated  perturba¬ 
tion  methods.  He  reported  on  the  use  of 
RNG  in  conjunction  with  a  large  eddy  sim¬ 
ulation  in  a  calculation  involving  32^ 
points  at  (wall  unit)  Reynolds  numbers 
of  460  to  1000.  (A  comparable  direct 
simulation  would  have  required  almost  10 
million  points!)  The  RNG  method  performs 
extremely  well  in  the  near-wall  region 
where  it  is  particularly  useful  in  un¬ 
steady  problems  where  the  log  law  cannot 
be  expected  to  hold.  In  such  situations 
the  RNG  theory  provides  excellent  predic¬ 
tion  of  such  quantities  as  the  phase  lag 
and  wall  shear  stress.  Excellent  predic¬ 
tions  of  the  convective  heat  transfer  at 
the  wall  can  also  be  obtained. 

Yet  another  approach  involves  mod¬ 
eling  turbulent  velocity  fields  with  sto¬ 
chastic  partial  differential  equatiorr. 
and  stochastic  boundary  conditions.  This 
approach  was  described  by  J.  Hunt  of  the 
University  of  Cambridge  (UK).  In  partic¬ 
ular,  Hunt  talked  about  the  so-called 
Lar.gevin  models  in  which  the  form  of 
the  particle  acceleration  spectrum  i«: 
provided  by  random  walk  or  random  Fou¬ 
rier  modes.  This,  Hunt  said,  provided 
a  "cheap"  turbulence  model,  and  he  de¬ 
scribed  conclusions  which  could  be  drawn 
from  such  calculations. 

6  I.ATTICE  GAS  METHODS 

Unlike  the  classical  approach  in 
which  the  fluid  is  treated  as  a  continu¬ 
um,  the  lattice  gas  method  represents  the 
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fluid  as  a  collection  of  discrete  fluid 
particles  moving  on  a  discrete  physical 
grid,  with  discrete  velocities  in  accord¬ 
ance  with  specific  interaction  rules  and 
conservation  laws.  In  other  words,  in¬ 
stead  of  writing  the  Navier-Stok.es  equa¬ 
tions  for  a  continuous  medium  and  then 
solving  them  by  a  discretization  process, 
one  obtains  discretized  equations  of 
motion  directly  by  discretizing  the  medi¬ 
um  and  its  associated  velocity  and  physi¬ 
cal  space  first. 

The  history  of  the  lattice  gas  meth¬ 
od  was  reviewed  by  J.-P.  Rivet  of  the  Ob- 
servatoire  de  Nice  (France).  In  the  first 
implementation  of  the  lattice  gas  method 
(Hardy  et  al.,  1973)  particles  were  con¬ 
strained  to  travel  on  a  rectangular  grid. 
This  was  a  natural  choice  given  the  ar¬ 
chitecture  of  the  specialized  computers 
(cellular  automata)  on  which  these  calcu¬ 
lations  were  carried  out.  However,  the 
resulting  stress  tensor  when  viewed  from 
a  macroscopic  perspective  failed  to  ex¬ 
hibit  the  isotropy  which  must  be  observed 
if  correspondence  with  the  Navier-Stokes 
equations  is  to  be  obtained.  In  1986, 
however,  Frisch  (Frisch  et  al.,  1986) 
discovered  that  the  desired  isotropy 
properties  could  be  recovered  if  a  tri¬ 
angular  (hexagonal)  instead  of  a  rectan¬ 
gular  grid  were  used  (see  Figure  2). 

Recent  theoretical  developments  have 
eliminated  other  restrictions  of  the 
early  lattice  gas  formulations.  Early 
models  were  unable  to  accommodate  thermal 


Figure  ?.  Interaction  Rules  of  the  KHP 
Model.  (Frisch,  et  al.,  1986) 


effects  (such  as  the  adiabatic  compres¬ 
sion  experienced  in  shock  waves)  due  to 
the  existence  of  a  single  particle  trans¬ 
lation  velocity  on  the  matrix.  The  devel¬ 
opment  of  more  flexible  collision  rules 
has  avoided  this  problem.  Another  ad¬ 
vance  has  been  the  calculation  of  three- 
dimensional  flows.  Up  until  the  past  few 
months,  this  seemed  to  he  impossible 
since  no  three-dimensional  lattice  was 
known  to  exhibit  the  required  isotropy. 
An  ingenious  solution  was  found  for  this 
problem  by  d'Humieres  (d'Humieres  et  al., 
1986)  in  which  a  three-dimensional  solu¬ 
tion  was  found  as  the  projection  of  an 
isotropy-preserving  four-dimensional  one. 

Another  recently  solved  problem  was 
that  the  lattice  gas  equivalent  of  the 
Navier-Stokes  equations  failed  to  exhibit 
Galilean  invariance.  That  is,  the  nature 
of  the  lattice  gas  equations  was  changed 
when  they  were  transformed  to  a  uniformly 
translating  frame  of  reference.  This  can 
be  traced  to  a  nonunitary  value  of  the 
parameter  which  is  used  to  scale  the 
values  of  time.  By  means  of  a  new  lattice 
gas  model  employing  "rest"  particles, 
d’Humi&res  and  Lallemand  (1986)  were  able 
to  increase  this  value  to  unity  and  thus 
restore  the  desired  invariance.  Of  a  more 
fundamental  natu re  appears  to  he  the  lim¬ 
itation  of  the  lattice  gas  method  to  low 
Mach  number,  low  Reynolds  number  flows. 
The  low  Mach  number  limitation  arises  due 
to  the  necessity  of  neglecting  terms 
which  depend  upon  the  Mach  number  if  cor¬ 
respondence  with  the  Navier-Stokes  equa¬ 
tions  is  to  be  achieved.  The  low  Reynolds 
number  limitation  results  from  the  neces¬ 
sity  of  limiting  the  number  of  particles 
participating  in  the  collision  processes 
in  order  to  have  a  workable  set  of  inter¬ 
action  rules.  Modest  increases  in  the 
Reynolds  number  are  possible  by  the  in¬ 
troduction  of  rest  particles,  however  at 
the  present  time  it  is  unclear  how  simu¬ 
lations  at  Mach  numbers  more  than  .»  few 
tenths  and  Reynolds  number  of  more  than  a 
few  thousand  can  be  achieved. 

Progress  in  the  past  ?  years  has 
been  spec  t  .icul  ar  with  lattice  gas  sim¬ 
ulations  lie  (ug  made  to  Investigate  buoy 
ancy  ellocts,  calculate  seismic  I’-waves, 
model  magnetohydrodynamic  flows,  solve 
reaction-diffusion  problems,  and  model 
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Figure  3.  I.attice  Gas  Simulations  of  the 
Flow  over  a  Backward  Facing  Step.  (Reyn¬ 
olds  Number  =  150) . 


interface  and  combustion  phenomena.  P. 
Lallemand  of  the  Laboratoire  de  Physique 
de  l'Ecole  Normale  (Paris)  described  sev¬ 
eral  recent  applications  including  the 
calculation  of  developing  Poisuille  flow 
at  the  entrance  of  a  rectangular  channel 
at  a  Reynolds  number  of  200  (see  Figure 
3) ,  flow  over  a  backward  facing  step  at 
a  Reynolds  number  of  150,  flow  around  a 
flat  plate  at  a  Reynolds  number  of  90, 
the  occurrence  of  Kelvin-Helmholtz  insta¬ 
bility  in  planar  shear  flow,  and  the  ap¬ 
pearance  of  Rayleigh-Tavlor  waves  in  a 
buoyancy-driven  flow. 


CONCLUSIONS 


The  contributions  which  applied 
mathematics  have  made  to  the  field  of 
computational  fluid  mechanics  particular¬ 
ly  in  the  fields  of  discrete  vortex  meth¬ 
ods,  multigrid  methods,  and  lattice  gas 
hydrodynamics  were  in  clear  evidence 
throughout  the  meeting.  Computational 
methods  have  become  a  very  popular  tech¬ 
nique  among  applied  mathematicians,  and 
numerous  fruitful  collaborations  between 
the  mathematical  and  fluid  mechanics  com¬ 
munities  have  resulted  in  the  development 
of  new  algorithms  and  the  development  of 
enhanced  convergence  procedures  for  ex¬ 
isting  ones.  Much  yet  remains  to  be  done, 
particularly  in  developing  innovative 
theories  (such  as  the  renormalization 
group  theory)  for  efficient  calculation 
of  turbulent  flows  and  methods  for  the 
calculation  of  turbulent  flows  with  re¬ 
action.  It  will  be  interesting  to  see 
how  much  progress  has  been  made  on  these 
topics  when  the  second  international  con¬ 
ference  on  industrial  and  applied  mathe¬ 
matics  is  held  in  4  years. 
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